Charge Transport Transitions and Scaling in Disordered Arrays of Metallic Dots 
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We examine the charge transport through disordered arrays of metallic dots using numerical simu- 
lations. We find power law scaling in the current-voltage curves for arrays containing no voids, while 
for void- filled arrays charge bottlenecks form and a single scaling is absent, in agreement with recent 
experiments. In the void-free case we also show that the scaling exponent depends on the effective 
dimensionality of the system. For increasing applied drives we find a transition from 2D disordered 
filamentary flow near threshold to a ID smectic flow which can be identified experimentally using 
characteristics in the transport curves and conduction noise. 
PACS numbers: 73.50.-h 



A wide variety of disordered systems exhibit threshold 
behavior and nonlinear response to an applied drive. Ex- 
amples include flux lines in disordered superconductors, 
[1-0] charge-density- waves (CDW's) pinned by impurities 
[I |, Wigner crystals in semiconductors with charge 
impurities, and colloids flowing over rough surfaces |Q. 
Another example is charge transport through metallic 
dot arrays. Middleton and Wingreen (MW) have con- 
sidered a model of this system in which the randomly 
charged dots are separated by tunnel barriers ||. They 
found threshold behaviors and scaling of the current- 
voltage curves of the form / = (V/Vr — 1) ■ In ID 
they obtain £ = 1.0, while for 2D they predict analyt- 
ically C — 5/3 and find in simulations £ = 2. For the 
2D systems the simulated flow patterns are not straight 
but form intricate meandering paths with considerable 
transverse fluctuations. These same types of meander- 
ing paths are also observed in the flow of flux-lines EL^J , 
Wigner crystals l|, and colloids JjJ above the depinning 
threshold. Experimental studies in metal dot arrays have 
also found scaling in the I-V curves for 2D and ID sys- 
tems |^-|l3|; however, the scaling exponents in these ex- 
periments exhibit a wide range of values. The studies for 
ID arrays find a scaling exponent of £ = 1.36, which 
is less than the value for 2D arrays predicted by MW, 
but still larger than the expected ID value of £ = 1.0. It 
is not known if this system is truly ID, or whether some 
meandering of the charge in 2D can still occur due to the 
finite width of the dots. It is also not known how the 
exponents would change (or whether there would even 
be scaling) upon changing the system from 2D to ID by 
gradually narrowing the array width. Other systems in 
2D exhibiting scaling near depinning also show a wide 
spectrum of scaling exponents (l]-||J|0, suggesting that 
the type of disorder and the effective dimensionality of 
the array plays a crucial role in the transport. 

Recently, to address the role of different types of dis- 
order on transport, Parthasarathy et al. ]l4| ] performed 
experiments on triangular monolayers of gold nanocrys- 
tals. Disorder is present in the form of charge disorder 
in the substrate as well as variations in the interparti- 



cle couplings. Structural disorder was also introduced 
by creating voids in the arrays. The void- free arrays ex- 
hibit robust power law scaling with £ = 2.25. A single 
power law could not be fit for the structurally disordered 
arrays. Parthasarathy et al. conjecture that in arrays 
with voids, the charge must be shuttled into bottleneck 
regions, reducing the amount of charge flow. In the theo- 
retical studies of MW where structural disorder was not 
considered, only power law behavior was observed 

Since the charge flow at depinning in the dot arrays re- 
sembles that seen in other systems such as flux lattices, it 
is of interest to ask whether some of the ideas developed 
in these other systems can be carried over to the metallic 
dot system. In disordered vortex systems, theory [p~5| p~7|| , 
simulation Jl5||l8| , and experiment |ll| show that at low 
drives the vortex flow is highly disordered, meanders in 
2D, and the overall lattice structure is destroyed. For 
high drives there can be a remarkable reordering transi- 
tion where the flow occurs in ID channels and the lattice 
regains considerable order. For strong quenched disorder 
in 2D, the highly driven phase forms a moving smectic 
with order in the transverse but not the longitudinal di- 
rection. Here the vortices move in well spaced channels, 
and the channels are decoupled from one another. 

In order to compare to the recent experiments and to 
explore the points raised above, we conduct simulations 
of a simple model for charge transport in 2D arrays with 
charge disorder and both with and without structural 
order. In our model we consider square and rectangular 
arrays of side N x M , with periodic boundary conditions 
in the x and, for the 2D systems, y directions, containing 
N c mobile charges. Under an applied drive fd a mobile 
charge on a site experiences a maximum threshold force 
fth before exiting the plaquette. For actual dot arrays 
the applied drive comes from an applied voltage V, and 
the energy to add an electron to a dot with charge q is 
Vth = <z/C where C is the capacitance of the dot. Charge 
flow will then occur for applied drives V > Vth for a 
single dot. In addition, the mobile charges interact via a 
Coulomb term U = q 2 jr. Since the Coulomb term is long 
range we use a summation technique |20| for numerical 
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FIG. 1. The scaling of the average velocity V vs applied 
drive f for (a) a 2D system with a fit of £ = 1.9 (solid line) 
and (b) a ID system with a £ = 1.0 fit (solid line). f t h is the 
threshold force at which depinning occurs, (c) A system of 
size 2 x 500, showing a fit with £ = 1.0. (d) A system of size 
5 x 500, with a £ = 1.41 fit. 

efficiency. For arrays without voids, we add disorder by 
selecting random thresholds Vth from a Gaussian distri- 
bution centered at V^ h . To study the effect of structural 
disorder, a portion P < 0.5 of randomly selected sites is 
effectively voided so that mobile charge cannot flow into 
them. We do not consider thermal effects since our sys- 
tem is in the Coulomb-blockade regime, so that charging 
energies are higher than the thermal energies. For in- 
creasing applied drive or voltage we measure the global 
charge flow or current /. We also measure the trajecto- 
ries of the flow and the fluctuations in the current, which 
is proportional to the conduction noise. 

We first consider the scaling for ordered arrays in 2D, 
ID, and finite width samples. In Fig. 1(a) we show the 
scaling of the average flow vs applied drive (I-V) curve 
for the 2D case, and in Fig. 1(b) the ID case. The curves 
in Fig. l(a,b) are for the largest systems we have consid- 
ered; for smaller systems the scaling regime is reduced. 
For 2D we find scaling with £ = 1.9 ± .15 in fair agree- 
ment with the simulation results £ = 2.0 of MW, but still 
lower than the £ = 2.25 found in Ref. [O for ordered ar- 
rays. The simulations for ID give a linear behavior for 
much of the curve, suggesting that if a scaling exponent 
could be ascribed it would be £ < 1.0; however, for drives 
near threshold, a single scaling cannot be applied and the 
curve bends up. We note that for the depinning of ID 
elastic objects such as CDW's one expects an exponent 
of C = 1/2 |Q. Experiments measuring the threshold of 
a single dot also find £ = 1/2 pj. Fig. 1(b) shows that 
the deviation from linearity occurs only very near thresh- 
old, so that the discrepancy between our results and MW 
may result from MW not being close enough to the 
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FIG. 2. The current paths for increasing applied drives 
in the 2D system shown in Fig. 1(a). Here (a) fd/ fth =1.1 
and (b) fd/ fth = 2.0 produce meandering 2D channels. In (c) 
fd/ fth = 10.0 produces straight ID channels. 

threshold to see the deviation. In Fig. 1(c) we show 
the scaling for a system with a finite width of 2 x 500, 
with £ = 1.0, and in Fig. 1(d) the scaling for a wider 
system of 5 x 500 is shown, where we find scaling with 
£ = 1.41. We note that experiments for the ID arrays 
Jiof find £ — 1.36. Our results suggest that scaling can 
occur for systems between 1 and 2 dimensions with the 
value of the exponent monotonically increasing from 1/2 
in ID to 2.0 in 2D. Additional evidence for the increase 
of the exponent as a function of dimensionality has been 
obtained in colbalt nanocrystal samples of finite thick- 
ness, with an effective dimensionality between 2 and 3. 
Here exponents of 2.2 < £ < 2.7 are observed [l2fl. Our 
results also suggest that the experiment in Jl4| is in an 
effective dimension higher than 2.0. All of the curves in 
Fig. 1 show a crossover to a linear regime at high drives, 
which was also seen in the earlier numerical work. We 
do not find any hysteresis for either the 2D or ID case. 
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FIG. 3. (a) Circles: dV/df of the curve shown in Fig. 1(a) 
for a 2D system. Squares: The noise power So as a function 
of drive, (b) Consecutive power spectra for increasing fd for, 
from bottom to top, f d /fth = 1-1, 2.0, 5.0, and 8.0. The 
curves have been shifted up for presentation. 

We next consider the flow transitions for increasing ap- 
plied drive in the 2D system. In Fig. 2 we show the cur- 
rent paths for three different applied drives above thresh- 
old. In Fig. 2(a), at fd/fth = 1-1 the flow follows me- 
andering paths in only a few regions of the sample, in 
agreement with the simulations of MW Q. These flow 
paths are stationary over time. For higher drives, as seen 
in Fig. 2(b) at fd/fth — 2.0, there is a crossover from 
the static paths to dynamic paths which open, close, and 
shift position over time. The filamentary flow in this 
case occurs everywhere in the system over time. For 
drives at and above the fd value where the I-V curve be- 
comes linear, as seen in Fig. 2(c) for fd/fth — 10.0, there 
is a crossover from the meandering 2D flow to straight 
ID ordered channels of flow. Fig. 2(c) shows that the 
charge moves only in ID channels without any jump- 
ing of charge between adjacent channels. The channels 
themselves carry different amounts of flowing charge due 
to the different average disorder along the rows. The 
charges in one channel do not synchronize with the flow 
of charge in adjacent channels; instead, the channels slide 
past one another. We term this a smectic flow state since 
the channels are periodically spaced in the transverse di- 
rection but are independently moving in the longitudinal 
direction. The transition from the disordered to partially 
ordered flow is very similar to the reordering transitions 
seen for driven vortex lattices 15 19| . We observe the 



same transitions in the systems of finite widths. 

In Fig. 3(a) we plot the dV/df curve, which is pro- 
portional to the resistance R, for the 2D system in 
Fig. 1(a). The crossover to the linear regime (with con- 
stant R) appears as the plateau region in dV/df. Also 
shown in Fig. 3(a) is the power So from one octave 
of the power spectra of the conduction noise at four 
different applied drives fd- So = J" a dvS{v), where 
S{v) = \JV x (t)e- i2nut dt\ 2 . The noise power shows a 
peak in the scaling regime near fd = 0.35 and then de- 
creases in the linear regime with a low value of Sq during 
the smectic flow. In Fig. 3(b) we show that a clearer sig- 
nature of the flow phases can be obtained by examining 
individual power spectra. For drives in the scaling regime 
we find a 1//" power spectra, which is indicative of the 
many different frequencies generated by the complex flow 
patterns illustrated in Fig. 2(a,b). The large noise power 
and 1// Q signals have also been associated with mean- 
dering disordered flow in superconducting vortices. For 
increasing drive a characteristic peak in the spectra be- 
gins to to appear and is most prominent in the smectic 
flow regime. This peak occurs when the charge in the 
smectic phase flows in ID paths along the dots which are 
in a periodic array of spacing a. The frequency at which 
the peak occurs is then v = v/a, where v is the average 
velocity of the charge. For perfectly ordered flow, the 
peaks would be very narrow. Since the channels have 
different amounts of flowing charge there is some disper- 
sion in the frequency. The presence of peaks in the power 
spectra suggests that it would be possible to observe an 
interference effect, or Shapiro steps, in the TV curves if 
an additional applied AC drive is imposed. The steps 
occur when the frequency of the AC drive matches the 
internal frequency of the system. 

Another experimental probe of the moving smectic 
phase is the presence of a transverse depinning barrier 
as first predicted in j[7) for elastic media. If a transverse 
force is applied to the already longitudinally moving sys- 
tem, then in the disordered regime there is no threshold 
for transverse motion and some charge will immediately 
begin to move in the transverse direction. For the high 
drive regime, after the ID moving channels have formed, 
there is a finite transverse threshold since the channels 
are effectively pinned in the transverse direction. There 
may also be interesting results for driving along different 
directions of the dot lattice. Although there is random- 
ness in the individual dot strength, the overall topological 
order of the array can break the symmetry so that certain 
directions may allow easier charge flow than others. 

We next consider the structurally disordered arrays. 
We find that for a fixed applied drive, the current is re- 
duced as the void fraction P is increased. This is un- 
derstandable since the charges must flow in increasingly 
winding patterns to pass the voids, which are effective 
obstacles. In Fig. 4(a) we plot V vs P up to P = 0.49 
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FIG. 4. (a) The average velocity at a fixed /<j = 0.2 for 
a system with increasing void fraction P. The solid line is a 
fit with V = 2(0.5 - P) 15 . (b) The velocity vs applied drive 
curve for a void fraction of P = 0.47. (c) The trajectories for 
a system with a void fraction of P = 0.47 at fd/fth = 3.0. 

for a fixed fd = 0.2. For the decreasing value of V, we 
find a best fit to a power law with V — Vb(0.5 — P) 1 - 5 . 
The current should go to zero for P = 0.5, or the per- 
colation limit, where for an infinite system voids would 
span the entire system. In our system there is still trans- 
port for P > 0.5 due to finite size effects. In Fig. 4(b) 
we show the current vs applied drive curve for a void 
fraction of P = 0.47. In this case, we cannot fit a single 
power law above threshold. In Ref. the additional 
features in the I-V curves for the structurally disordered 
arrays were conjectured to occur due to bottleneck effects 
caused by the void regions. In Fig. 4(c) we illustrate the 
current paths for fd/fth — 3.0 and P = 0.47, showing 
that there is considerable flow through the system, but 
that in certain well-defined areas no flow is occurring. 
Averaging the trajectories over a longer time produces 
the same flow patterns shown in Fig. 4(c). This is in 
contrast to the flow pattern in Fig. 2(b), which changes 
over time, so that for long times flow occurs in all regions 
of the sample. In Fig. 4(c) some bottlenecks can also be 
seen in the form of regions where the trajectories are 
compressed. For the structurally disordered arrays the 
transition to the smectic flow at higher drives is absent 
since straight ID flows cannot occur even in the ohmic 
regime. The conduction noise shows the same l/f a be- 
havior as the ordered arrays but the peaks in the noise 
spectra are absent in the ohmic regime. 

In summary, we have investigated charge transport in 



structurally ordered and disordered arrays. For ordered 
arrays we find scaling in the current vs applied drive 
curves with ( — 1.9 in 2D. For ID arrays we find that the 
scaling exponent is less than one. Scaling still occurs for 
systems with finite width, with the exponent increasing 
toward £ = 2.0 for increasing sample widths. For increas- 
ing applied drive in 2D, we show that the crossover to 
ohmic behavior coincides with a change in the flow from 
2D meandering to straight ID channels or smectic flow. 
Evidence for this change in the flow also appears in the 
form of a crossover in the power spectra, which shows a 
broad 1// Q signature in the disordered flow regime, and 
a characteristic peak or washboard signal in the smectic 
flow regime. For disordered arrays where a fraction of the 
sites are replaced with voids, a single power law cannot 
be fit to the I-V curve in agreement with recent experi- 
ments. The transition from the 2D disordered flow to the 
ID channel flow is absent in the structurally disordered 
arrays. 
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